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Abstract. Observations of kHz quasi-periodic oscillations (QPOs) in the X-ray fluxes of low-mass X-ray binaries 
(LMXBs) have been used in attempts to constrain the external metric of the compact members of these binaries, 
as well as their masses and the equation of state of matter at supranuclear denisties. We compute the maximum 
orbital frequency of stable circular motion around uniformly rotating strange stars described by the MIT bag 
model. The calculations are performed for both normal and supramassive constant baryon mass sequences of 
strange stars rotating at all possible rates. We find the lower limits on the maximum orbital frequency and discuss 
them for a range of masses and for all rotational frequencies allowed in the model considered. We show that for 
slowly and moderately rotating strange stars the maximum value of orbital frequency can be a good indicator of 
the mass of the compact object. However, for rapidly rotating strange stars the same value of orbital frequency 
in the innermost stable circular orbit is obtained for stars with masses ranging from that of a planetoid to about 
three solar masses. At sufficiently high rotation rates of the strange star, the rotational period alone constrains 
the stellar mass to a surprisingly narrow range. 
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1. Introduction 

The recently discovered kHz quasi periodic oscillations 
may be used as a probe of the inner regions of accretion 
disks in compact stars and hence also of some of the prop- 
erties of the central object. For example, these oscillations 
have been used to derive estimates of the mass of the cen- 
tral source in some sources (e.g., Kaaret et al., 1997; Zhang 
et al., 1998; Kluzniak 1998). All these estimates assume 
that the maximum observed QPO frequency is the orbital 
frequency at the innermost stable circular orbit (ISCO), 
as suggested by Kluzniak et al. (1990). 

In the first study of the possible presence of strange 
stars in those LMXBs which exhibit kHz QPOs, Bulik et 
al. (1999a, b) showed that slowly rotating strange stars, 
described by the simple MIT bag model with massless 
and non-interacting quarks, have orbital frequencies at 
the marginally stable orbit higher than the maximum fre- 
quency of 1.07 kHz in 4U 1820-30, reported by Zhang et 
al. (1998). However, the ISCO frequencies can be as low 
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as 1 kHz when more sophisticated models of quark mat- 
ter (MIT bag with massive strange quarks and lowest or- 
der QCD interactions) and/or rapid stellar rotation are 
taken into account (Stergioulas et al., 1999; Zdunik et al., 
2000a, b). The main result of these considerations was that 
the lowest orbital frequency at the ISCO is attained either 
for non-rotating massive configurations close to their max- 
imum mass limit, or for configurations at the Keplerian 
mass-shedding limit. In order to obtain an orbital fre- 
quency in the ISCO as low as 1.07 kHz for slowly rotating 
strange stars one has to consider a specific set of (rather 
extreme) MIT bag model parameters (Zdunik et al. 2000a, 
b) allowing gravitational masses higher than 2.2 Mq. In 
contrast, for stars rotating close to the Keplerian limit, a 
low orbital frequency at the ISCO can be obtained for a 
broad range of stellar masses (Stergioulas et al., 1999). 

A new description of quark matter (Dey et al., 1998) 
has been used for constructing numerical models of static 
strange stars and strange stars rotating with two different 
fixed frequencies 200 Hz and 580 Hz, and for calculating 
frequencies of marginally stable orbits around these stars 
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(Datta et al., 2000). The authors conclude that very high 
QPO frequencies in the range of 1.9 to 3.1 kHz would im- 
ply the existence of a non-magnetized strange star in an 
X-ray binary rather than a neutron star. Gondek-Rosinska 
et al. (2001) calculated the maximum orbital frequencies 
for the Dey et al. quark stars rotating at various rota- 
tion rates and showed that the rotational and maximum 
orbital frequencies for these stars are indeed much higher 
than those for neutron star models and for strange stars in 
the MIT bag model. The maximum orbital frequencies for 
the compact objects modeled with the Dey et al. (1998) 
equation of state are > 1.5 kHz, and always higher than 
the kHz QPO frequencies observed to date. 

Recently, new complications in the problem were 
pointed out. The innermost stable circular orbit exists also 
in Newtonian gravity (Kluzniak, Bulik, Gondek-Rosihska 
2000; Amsterdamski et al., 2000; Zdunik & Gourgoulhon 
2001). The presence of a crust has been investigated 
for normal evolutionary sequences by Zdunik, Haensel 
& Gourgoulhon (2001). Typically the crust increases the 
maximum orbital frequency at the Keplerian limit. 

In this paper we consider constant baryon mass se- 
quences of strange stars described by the simplest MIT 
bag model (with massless and non-interacting particles) 
for various values of the baryon mass and calculate the 
frequency in the co-rotating ISCO (if it exists) for each se- 
quence. We find and discuss lower limits to the maximum 
orbital frequency (taking into account various secular and 
dynamical stability limits, as well as physical constrains 
relevant for accretion in LMXBs). We then extend our re- 
sults to all MIT equations of state. Our calculations are 
described in section 2 and the results are presented in 
section 3. Finally, section 4 contains a summary of our 
results. 

2. The model star calculations 

Typically, strange stars are modeled (Alcock et al., 1986; 
Haensel et al., 1986) with an equation of state based on 
the MIT-bag model of quark matter, in which quark con- 
finement is described by an energy term proportional to 
the volume (Farhi and Jaffe 1984). The equation of state 
is given by the simple formula 

P = a(p-po)c^ (1) 

where P is the pressure, p the mass-energy density, and 
c the speed of light. It was shown by Zdunik (2000) that 
all MIT bag models can be approximated very well by 
equation (1), where a is a constant of model-dependent 
value in the range 0.289 (for a strange quark mass of ms= 
250 MeV) to 1/3 (for massless quarks). In the numeri- 
cal calculations reported in the present paper we describe 
strange quark matter using the simplest MIT bag model 
(with massless and non-interacting quarks), accordingly 
we take a = 1/3 and po = 4.2785 x 10^^ g cm~^ = pad 
(corresponding to a bag constant of S = 60 McV fm~^). 
In general, eq. (1) corresponds to self-bound matter with 



density po at zero pressure and with a fixed sound ve- 
locity (i/ac) at all pressures. For a fixed value of a, all 
stellar parameters are subject to scaling relations with 
appropriate powers of po (see, e.g., Witten 1984; Zdunik 

2000) . Approximate scalings with a were considered by 
Stergioulas et al. (1999) for the limiting case of rotation 
at Keplerian frequencies. Although a precise scaling-law 
with a for other cases has not been derived, we will show 
that the lower limits on the maximum ISCO orbital fre- 
quency obtained for the simplest MIT bag model can eas- 
ily be used to find lower limits on the maximum orbital 
frequency for any other value of a and po- 

We have calculated the maximum orbital frequency 
around uniformly rotating strange stars using two differ- 
ent highly accurate relativistic c;odes. The first code is 
based on the multi-domain spectral methods developed by 
Bonazzola et al. (1998). This method has been used previ- 
ously for calculating properties of rapidly rotating strange 
stars described by the MIT bag model (Gourgoulhon et 
al., 1999) and by the Dey model (Gondek-Rosinska et al., 
2000a). This code was also used for calculating the max- 
imum orbital frequency around normal sequences of MIT 
strange stars (Zdunik et al., 2000b, Zdunik & Gourgoulhon 

2001) and of the Dey model strange stars (Gondek- 
Rosinska et al., 2001). The second code (Stergioulas & 
Friedman 1995, see Stergioulas 1998 for a description) has 
been used to compute accurate constant baryon mass se- 
quences and models at the Keplerian mass-shedding limit, 
where the first code has difficulties converging to a solu- 
tion. In this code the equilibrium models arc obtained fol- 
lowing the KEH method (Komatsu et al., 1989), in which 
the field equations are converted to integral equations us- 
ing appropriate Green's functions. Models computed by 
both codes agree very well in all computed properties. 

We calculate the maximum frequency of a co-rotating 
test particle in a stable circular orbit in the equatorial 
plane about a strange star, /orb, max. By testing stability 
of the orbital motion, we determine whether or not stable 
orbits extend to the surface of the star (at radius i?*). If 
they do. the maximum orbital frequency is the keplerian 
frequency at the surface, /orb, max = /orb(-R*)- If they do 
not, we find the radius of the innermost marginally stable 
orbit, i?ms, and its frequency /ms (see, e.g. Friedman, Ipser 
& Parker 1986). In this latter case, the orbital frequency 
in the ISCO is the maximum circular orbital frequency, 
/orb, max = /ms • This frequency can then be compared 
with the observed kHz QPOs, having in mind some model 
of QPOs relating their frequency to the orbital frequency 
and the location of the ISCO (e.g., Kluzniak et al., 1990, 
Abramowicz and Kluzniak 2001, Wagoner et al., 2001). In 
so doing, we neglect the effects of fluid pressure, magnetic 
fields, radiation drag, etc., on the location of the innermost 
circular orbit in the accretion disk. 

We construct equilibrium sequences of rotating com- 
pact strange stars with constant baryon mass. We com- 
pute both normal and supramassive sequences. A sequence 
is called normal if it terminates at the zero angular mo- 
mentum limit with a static, spherically symmetric solu- 
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Fig. 1. The dependence of the maximum orbital frequency on the frequency of rotation ri/27r. Each thin solid line 
corresponds to a sequence of equilibrium models with a constant baryon mass, and is labeled by this baryon mass 
in solar mass units, as well as (in parentheses) by the gravitational mass of the static configuration, if it exists. The 
dashed line indicates stars marginally stable to axisymmetric perturbations. The angular momentum increases along 
each sequence from J = for static configurations, or Jmin on the dashed line for supramassive stars, to Jmax for the 
Keplerian configurations (thick solid line). One sequence for a very low mass star (with Mb = O.OIMq) is shown, the 
critical point on this sequence for Newtonian dynamical instability to non-axisymmetric perturbations is indicated by 
an asterisk and dynamically unstable configurations are denoted with the dotted line. The hatched area corresponds to 
a region excluded for rotating strange stars in the model considered. The inset shows in detail how the supramassive 
sequence reaches the stability limit to quasi-radial oscillation. 



tion, and it is called supramassive if it terminates at the 
axisymmetric stability limit, instead. The boundary be- 
tween these two sequences is the sequence with the max- 
imum baryon mass of a static configuration. The angular 



momentum (and the central density) of a star changes 
monotonically along each sequence. 

For a supramassive constant baryon mass sequence 
a model is stable to axisymmetric perturbations, if 
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f dJ\ ^ , fdM\ , , . , 

— — < (or — — < 0) where J is the an- 

gular momentum, M is the gravitational mass, Mb is the 
baryon mass of the star and pc is the central density of 
the star (Friedman, Ipser & Sorkin 1988; Cook, Shapiro 
& Teukolsky 1994). 

The star in a sequence reaches the Keplerian limit (the 
mass-shedding limit) if the velocity at the equator of a 
rotating star is equal to the velocity of an orbiting particle 
at the surface of the star. 



3. Results 

The properties of the ISCO around strange stars differ 
from those around neutron stars (Miller, Lamb & Cook 
1998) — for strange (i.e., quark) stars of any mass the 
ISCO always exists at sufficiently rapid stellar rotation 
rate (Stergioulas et al., 1999; Zdunik et al., 2000b), and 
for strange stars of moderate and high baryon masses the 
ISCO always exists for any rotation rate (Zdunik et al., 
2000b; Gondck-Rosinska et al., 2001). For lower mass stars 
(e.g., the I.QMq baryon mass sequence in Fig. 1) stable 
orbits extend down to the stellar surface for moderate ro- 
tation rates, but the gap is present for both slowly and 
rapidly rotating stars. For lower masses yet, the gap is 
present only at high rotation rates and the ISCO is present 
even if relativistic effects are negligible (as in substellar- 
mass configurations, e.g., for M < O.OlAfo). In the 
Newtonian limit, the gap between the marginally stable 
orbit and the stellar surface is produced by the oblateness 
of the rapidly rotating low-mass quark star (Kluzniak, 
Bulik & Gondek-Rosinska, 2000; Amsterdamski et al., 
2000; Zdunik & Gourgoulhon 2001). For Keplerian config- 
urations the ISCO is always well above the stellar surface 
for all possible masses. 

We present the dependence of the maximum orbital 
frequency on the frequency of rotation /rot = fi/27r for 
a set of constant baryon mass sequences of strange stars 
in Figure 1. In Figure 2 we show how the gravitational 
mass depends on rotational frequency for the consid- 
ered sequences. The normal sequences start at /rot = 0, 
and continue with increasing angular momentum up to 
the Keplerian configuration. As it turns out, the very 
low mass limit, represented in Figure 1 by the case of 
Mb = O.OIMq, is described well by a uniform-density 
fluid, i.e., the Maclaurin spheroid (Chandrasekhar 1969). 
For this sequence, an ISCO exists only for rotation rates 
larger than ^ 970 Hz, and we indicate the onset of the 
(Newtonian) dynamical instability to non-axisymmetric 
perturbations by an asterisk. 

The intermediate-mass case is shown by the constant 
baryon mass sequence with the baryon mass Mb = I.GMq, 
for which the ISCO does not exist between the rotation 
rates with frequency of /rot ~ 300 Hz and /rot ~ 800 Hz 
(see also Zdunik et al., 2000b). 

More massive stars in normal sequences are shown for 
the cases of Afb = 2.OAf0 and 2.3M0. In these cases, 
the marginally stable orbit exists along the entire con- 




Fig. 2. Gravitational mass as a function of the rotational 

frequency for the same constant baryon mass sequences as 
shown in Figure 1. The symbols are the same as in Fig 1. 
The short dashed line corresponds to the upper limit on 
gravitational mass for the model considered. 



stant baryon mass sequence. The ISCO frequency ini- 
tially increases with rotational frequency, as calculated in 
the slow-rotation approximation by Kluzniak & Wagoner 
1985. When the small-rotation-rate assumption is no 
longer satisfied, the ISCO frequency reaches a maximum 
and begins to decrease mainly due to the rotational defor- 
mation of the star (Sibgatulin & Sunyaev 1998; Shibata 
& Sasaki 1999). The sequence terminates when the star 
has reached the Keplerian configuration (mass-shedding 
limit). An extreme case of normal sequences is the se- 
quence starting with the maximum mass static configura- 
tion of Mb = 2.624Mq (indicated as 2.6 M© in Figures 1 
and 2). 

Supramassive sequences start at the onset of the in- 
stability to axisymmetric perturbations (thick dashed line 
in Figure 1 and in subsequent figures.) Sequences with 
Mb = 2.7 Mq and 2.8M0 are examples of low mass supra- 
massive sequences. The behavior of the maximum orbital 
frequency vs. frequency of rotation is basically similar to 
the normal sequences described above. 

Finally, we present the sequences of models with Mb = 
3.2, 3.4, and 3.6M0. For these models of extremely high 
mass, the rotational frequency decreases with increasing 
angular momentum because the stars are initially highly 
deformed. The ISCO frequency only decreases along these 
sequences until the stars reach the mass shedding limit. 
The configuration with the highest rotation frequency is 
not Keplerian. For a high-mass supramassive strange star 
sequence (e.g., Mb = 3.2), the Keplerian configuration is 
the one with the lowest rotational frequency in the se- 
quence (see Fig. 2). 

Note that the rotational frequency does not necessar- 
ily change monotonically along a constant baryon mass 
sequence, since a star changes its shape and, consequently. 
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Fig. 3. The ISCO frequency as a function of central den- 
sity for sequences with a fixed rotational frequency, as 
labeled. The short dashed line denotes the configurations 
with the maximum possible gravitational mass for given 
frequency (see Figure 2), the dash-dotted line corresponds 
to stars with maximum baryon mass, and the thick dashed 
line represents the axisymmetric stability limit. 

its moment of inertia with increasing angular momentum 
(see also Cook, Shapiro & Teukolsky 1992; Zdunik et al., 
2000b; Gondek-Rosiiiska et al., 2GG0a). For a lower-mass 
supramassive sequence and a normal sequence the rota- 
tional frequency decreases with increasing angular mo- 
mentum just below the mass-shedding limit (character- 
istic turning back in Figure 1 and 2 for each sequence). 
A similar behavior is observed in the case of models with 
small mass, e.g. for the constant baryon mass sequence 
with Mb = O.O1M0. Thus, the effect of spin-down with 
increasing angular momentum is of Newtonian origin, and 
is related to the change of shape of the star. 

In Figure 1 one can see that the value of the ISCO fre- 
quency alone is not a sufficient indicator for the mass of 
the compact object. For a rapidly rotating strange star the 
same value of orbital frequency at the innermost stable cir- 
cular orbit, e.g. 1.2 kHz, can correspond to a range of stars 
differing in mass by a factor of a hundred, or even more, 
as stressed already by Kluzniak, Bulik & Gondek-Rosihska 
2000. However, it is worth noting that for rotational fre- 
quencies exceeding about 1.3 kHz, i.e., for rapidly rotating 
supramassive stars, the rotational period alone constrains 
the stellar mass to a narrow range. 

3.1. Lower limits on the ISCO frequency 

Figure 1 shows that at any given stellar rotational fre- 
quency the maximum orbital frequency /orb,max(/) is hm- 
ited from below by the least frequency among the follow- 
ing three "extreme" sequences of stars: a) the marginally 
stable configuration to axisymmetric perturbations (thick 
dashed line); b) configurations rotating at Keplerian fre- 



quencies (thick solid line) and c) the Newtonian limit 
of low-mass stars, tending to Maclaurin spheroids when 
A/ ^ 0, J — > 0, p ^ Pq (solid line ending with the aster- 
isk). For all these limiting cases, an ISCO exists. The max- 
imum orbital frequencies are also bounded from above by 
a fourth limiting line (not shown) , which will be discussed 
elsewhere: the envelope of highest values of /orb. max- 

a) We find that at a given low or moderate value of stel- 
lar rotational frequency the lowest ISCO frequency is ob- 
tained for the configurations at the axisymmetric instabil- 
ity limit. As was shown by Gourgoulhon et al. (1999) and 
Gondek-Rosihska et al. (2000a) a supramassive strange 
star, just as a supramassive neutron star (Cook et al., 
1992), prior to its collapse to a black hole increases its 
spin as it loses angular momentum. In such cases, two 
models with the same baryon mass but with different cen- 
tral density and angular momentum can have the same 
rotational frequency (see inset in Figure 1). Along a se- 
quence of fixed rotational frequency, the model with the 
lowest ISCO frequency is not the model with the high- 
est gravitational or highest baryon mass. The dependence 
of the ISCO frequency as function of the stellar central 
density for sequences of stars with a fixed rotational fre- 
quency is shown in Figure ^. As the central density in- 
creases, the ISCO frequency decreases and the sequence 
passes through the configuration with maximum gravita- 
tional mass (denoted with short dashed line in Figure |^ 
and Figure 2), and then the maximum baryon mass config- 
uration (denoted with the dash-dotted line) . Before reach- 
ing the minimal value of /orb, max, the sequence reaches 
the limit of stability to axisymmetric perturbations (thick 
dashed line). Thus, for a fixed rotational frequency, the 
lowest possible maximum orbital frequency is attained for 
the model which is marginally stable to gravitational col- 
lapse. 

b) At high rotation rates the lowest ISCO frequency 
is attained for the Keplerian (mass-shedding) configura- 
tions (thick solid line in Figure 1). The stellar configu- 
rations at this limit are highly deformed. We calculated 
Keplerian configurations beginning with the fastest ro- 
tating stable configurations (with Mb = 3.76Mq), end- 
ing with a strange star with Afb = 0.4Mo. Note that 
for rapidly rotating (close to the Keplerian limit) strange 
stars, the ISCO frequency is smaller than the rotational 
frequency. 

c) In the small-mass limit the ISCO only exists at high 
stellar rotation rates. This was discovered numerically, and 
the Newtonian origin of the ISCO in this case was first un- 
derstood by comparison with toy models of disk-like dis- 
tributions of matter (Kluzniak, W., Bulik, T.,& Gondek- 
Rosihska 2000; Zdunik & Gourgoulhon 2001). In the 
Newtonian limit {M — > 0, J/M'^ oo) strange stars tend 
to a nearly uniform-density configuration {p po), and 
these can be described well by the Maclaurin spheroids, for 
which analytic formulae for orbital frequencies have been 
derived and found to agree with the numerical models 
(Amsterdamski et al., 2000). The functional dependence 
of the maximum orbital frequency on the stellar rotation 
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Fig. 4. The effect of non-zero strange quark mass on 
the lower limits on the maximum orbital frequency. The 
dependence on po has been scaled out. Here, the thin 
dotted-dashed lines correspond to the limits for the case of 
m.s = 250 MeV (a=0.289), and all other lines and symbols 
have the same meaning as in Figure 1, which was obtained 
for a = 1/3 ( uis = 0). 

rate is nearly universal for all strange stars with masses 
lower than 0.1 Mq, in the model considered here. 

3.2. The lower limits on ISCO for other MIT bag 
models 

Figure 1 and the discussion above assumed the simplified 
MIT bag model with massless and non- interacting quarks. 
The results can easily be scaled for any value of the density 
at zero pressure, with the use of the scaling / cx pj^^, 
where / denotes either of the frequencies (rotational or 

— 1/2 

orbital), and M cx . We can easily find lower limits 

on ISCO for any other MIT bag model since Zdunik (2000) 
showed that all MIT bag model equations of state can be 
well approximated by eq. (1). 

In order to find the dependence of the lower limits on 
a we display in Figure || the limits calculated for mg = 
and for uis = 250 MeV. We have chosen these values of the 
strange quark mass because it corresponds to the highest 
value a = 0.333, and the lowest value of a = 0.289 for 
the MIT bag model equation of state. The dependence 
on the density of quark matter at zero pressure has been 
scaled out by multiplying the frequencies with the square 
root of pfid/po, where our fiducial value of po was used: 
pfid = 4.2785 X lO^'' g cm~^. The numerical values of the 
frequencies in Fig. 4 correspond to the value of the bag 
constant B = 60 MeV fm~^ used throughout this paper. 

The limit on the maximum orbital frequency, set by 
the low mass stars does not depend on the value of a. 
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Fig. 5. The ratio of the kinetic energy to the absolute 
value of the gravitational potential energy (defined as in 
Friedman et al., 1986) as a function of the rotational fre- 
quency for several constant baryon mass sequences of MIT 
bag model strange stars. The meaning of lines is the same 
as in Fig 1. The critical value of for the onset of 

dynamical stability in the low mass limit is the Newtonian 
value of 0.27 indicated by an asterisk. 

The limit at rapid rotation rates depends only slightly 
on the value of a. The main difference appears in the 
slow and moderate rotation regions, where the lower 
limit is set by the sequence of marginally stable stars. 
However, here one can use the fact that the limit is only 
shifted along the vertical axis. We find that the formula 
Jorb,max (flj/rot) = 0.3 /rot + /orb,max(a, 0) is a good ap- 
proximation for describing the lower limit on the ISCO 
frequency for slowly and moderately rotating strange stars 
described by any MIT equation of state. Here, the maxi- 
mum orbital frequency of the nonrotating model is given 
by its Schwarzschild-metric value 

/orb,max(a,0) = 2.2 kHz/(A/(a)/MQ) 

(e.g., Kluzniak et al., 1990), and M(a) is the maximum 
mass in general relativity of a static configuration of the 
(quark) fluid described by eq. (1) — in short, the maximum 
mass of a non-rotating strange star. 

3.3. Factors that could alter the lower limits 

There are several factors which could affect the lower lim- 
its considered in the previous section. We discuss here: i) 
different kinds of instabilities that could operate in rapidly 
rotating strange stars, ii) the existence of a solid crust on 
the accreting strange stars, iii) details of the accretion- 
induced spin-up of strange stars in LMXBs. 

i) Rapidly rotating compact stars may be subject to 
non-axisymmetric rotational instabilities. A widely used 
indicator for the onset of instability is the ratio of the 
kinetic energy to the absolute value of the gravitational 
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potential energy T/\W\. We plot the value of T/\W\ (de- 
fined as in Friedman et al., 1986) for several constant 
baryon mass sequences in Figure |[ The values of r/|Vl^| 
for strange stars rotating at the Keplerian velocity are 
very large and can be even higher than = 0.32.|^ 

The is very high for all linear self-bound equations 

of state (see Gondek-Rosinska et al., 2000b for different 
MIT bag models and Gondek-Rosinska et al., 2000a for 
the strange stars described by the Dey model). The val- 
ues of r/|M^| for strange stars are much larger than for 
any neutron star models (for which the maximum T/|VF| 
ranges typically from 0.1 to 0.14, (see e.g. Nozawa et al., 
1998). In contrast to neutron stars, the T/|Vl^| for strange 
stars increases with decreasing stellar mass for all models. 
Very large values of T/|Vl^| are characteristic not only for 
supramassive stars but also for rapidly rotating strange 
stars with small and moderate masses. The large value of 
T/|W| results from a flat density profile combined with 
strong equatorial flattening of rapidly rotating strange 
stars. No stability analysis in general relativity has been 
performed yet for strange stars. 

There exists two different classes and corresponding 
time-scales for non-axisymmetric instabilities: a dynami- 
cal instability, growing on the hydrodynamical time-scale 
and secular instabilities (due to dissipative mechanisms) 
which grow on much longer time-scales. In Newtonian the- 
ory, rapidly rotating Maclaurin spheroids are secularly un- 
stable to a bar-mode formation, if T/|l^| > (T/|Vt^|)crit = 
0.1375 (Jacobi/Dedekind bifurcation point) and dynam- 
ically unstable if {T/\W\ > 0.2738). A secular instabil- 
ity can grow only in the presence of dissipative mech- 
anisms, such as shear viscosity (Roberts & Stewartson 
1963) or gravitational radiation (the Chandrasekhar- 
Friedman-Schutz instability (hereafter CFS instability; 
Chandrasekhar 1970; Friedman & Schutz 1978; Friedman 
1978). If both viscosity and gravitational radiation are op- 
erative, they act against each other. Which of the two 
dissipative mechanisms dominates is highly dependent on 
the temperature of the star. By a suitable choice of the 
strength of viscosity relative to gravitational radiation, 
it is possible to stabilize the Maclaurin sequence up to 
the point of dynamical instability (Shapiro &; Teukolsky 
1983). 

In general, the viscosity-driven instability could be op- 
erating in accreting (old and cold) compact objects, while 
the CFS instability could operate in both newly-born com- 
pact stars and in LMXBs. The most important pulsation 
modes for the CFS instability are the r- and /-modes (see 
Stergioulas 1998, Andersson & Kokkotas 2001, Friedman 
and Lockitch 2001, for recent reviews of these instabili- 
ties). 

The dynamical instability against bar-mode deforma- 
tion of rapidly rotating polytropic stars was studied by 

^ In Gourgoulhon et. al., 1999 (see their Fig. 7), T/\W\ < 
0.24, because configurations denoted as Keplerian, were in fact 
configurations at which the numerical procedure stopped con- 
verging. 



Shibata, Baumgarte & Shapiro (2000) in full general rela- 
tivity and by Saijo et al., 2001 in the first post- Newtonian 
approximation. The authors considered differcntionally 
rotating proto-neutron stars. Their main result is that in 
general relativity the onset of dynamical instability oc- 
curs for a somewhat smaller critical value of (T/|VF|)crit ~ 
0.24 — 0.25 than the Newtonian value ~ 0.27 for incom- 
pressible Maclaurin spheroids. 

The few numerical studies of secular instabilities in 
rapidly rotating stars, in general relativity show that non- 
axisymmetric modes driven unstable by viscosity no longer 
coincide with those driven unstable by gravitational ra- 
diation. The (T/|VF|)crit depends not only on the dissi- 
pative mechanism but also on the compaction parameter 
M / i? of a star. The viscosity-driven instability for neutron 
stars has been studied in general relativity by Bonazzola, 
Frieben & Gourgoulhon (1996, 1998), who find that it 
sets in for r/|II^| somewhat larger than the Newtonian 
value of (T/|VF|)crit = 0.1375. A stronger weakening of 
the viscosity-driven instability by general relativity is sug- 
gested by a post-Newtonian study of the viscosity-driven 
instability in incompressible stars by Shapiro & Zane 
(1998), who find that the (T/|PF|)crit could be as large as 
0.25 for M/R ^ 0.2. On the other hand, the gravitational- 
radiation driven instability is strengthened by relativity 
and the bar {I = 2) /-mode can become unstable for 
typical neutron star equations of state (Stergioulas & 
Friedman 1998). For a wide range of realistic equations of 
state, Morsink, Stergioulas & Blattning (1999) find that 
the critical for the onset of the instability in the 

I = 2 /-mode is only ~ 0.08 for models with mass 

1.4M0 and T/|W^| ~ 0.06 for maximum mass models. 

The CFS-instability of r-modes in strange stars has 
been studied by Madsen (2000), in the Newtonian limit. 
According to this study, the instability can operate at a 
minimum rotational frequency of ~ 400Hz and at a tem- 
perature of T - lO^K. Since strange stars can, at most, 
have only a thin solid crust, they cannot store heat in 
the interior during accretion. Thus, accreting strange stars 
are likely to have interior temperatures comparable to 
T ~ lO^K. If the r-mode instability is not damped by 
other mechanisms then the above results seem to rule out 
the presence of low and moderate mass stars with a nor- 
mal quark matter in LMXBs (if the observed 1.07 kHz 
QPO is the orbital frequency at the ISCO). 

If, contrary to our current understanding, the r-mode 
instability is not effective in accreting strange stars, then 
the /-mode instability or the viscosity-driven instability 
could be more relevant. More detailed studies of these 
instabilities specifically for relativistic strange stars are 
needed. The dynamical instability, even if present at 
r/|W^| ~ 0.24, only marginally affects the limits presented 
in Fig. 1. 

ii) Another factor affecting the lower limits presented 
in Figure 1 would be the presence of a solid crust. 
This effect has been discussed by Zdunik, Haensel & 
Gourgoulhon (2001). The presence of the crust will not 
affect the lower limit in the small and moderate rotation 
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case. Here the stability of the strange stars is basically 
determined by the stability of the strange core and the 
properties of the crust are negligible. However, for rapid 
rotation, when the lower limit on the maximum orbital 
frequency is determined by the Keplerian configurations, 
the presence of the crust will increase the lower limit. The 
crust is of low density and for rapidly rotating strange 
stars it becomes extended in the equatorial region. Thus 
the mass-shedding limit is reached at lower rotation rates 
than in the case of bare strange stars. Zdunik, Haensel 
& Gourgoulhon (2001) showed that the values of max- 
imum orbital frequency for strange stars with maximum 
crust (up to the neutron drip point), rotating at the mass- 
shedding limit, arc about 20% larger compared to case of 
bare strange stars. Also, for strange stars with a crust, the 
ISCO exists only for models with large mass. 

iii) When a strange star is spun-up by accretion of 
material that leaves the ISCO and plunges onto the sur- 
face of the star, the largest attainable rotation rate is 
not precisely the rotation rate at the mass-shedding limit. 
Material that leaves from the ISCO arrives at the sur- 
face of a strange star after following a path in which it 
conserves angular momentum. At the surface, the mate- 
rial has angular velocity larger than what it had at the 
ISCO, but still somewhat smaller than the angular veloc- 
ity needed to spin-up the star to the angular velocity of 
the equilibrium configuration at the mass-shedding limit. 
The difference between the maximum rotation and mass- 
shedding configurations is, however, small. 

4. Summary 

Relativistic orbital frequencies are crucially related to the 
spin-up of neutron stars and black holes and to high en- 
ergy accretion phenomena in low mass X-ray binaries, in- 
cluding quasi-periodic oscillations in the X-ray flux. We 
compute the maximum orbital frequency of stable circular 
motion around strange stars, and present it as a function 
of the stellar rotation rate for normal and supramassive 
constant baryon mass sequences of stars rotating at all 
possible rates. We find that in contrast to neutron stars, 
the mass of a rapidly rotating strange star cannot be even 
approximately inferred from the orbital frequencies alone. 
For strange stars described by the MIT bag model, the 
same frequency in the innermost stable orbit (e.g., 1.2 
kHz) is obtained for strange stars with rotational periods 
ranging from infinity to about 0.6 ms, and the mass rang- 
ing from that of a planetoid to about three solar masses. 
One should note however that the rapidly rotating stars 
in this region are quite likely to be unstable. 

We find that for slow to moderate rotation rates the 
lower limit on the maximum orbital frequency around 
strange stars is attained for configurations marginally sta- 
ble to axisymmetric perturbations, while at rapid rotation 
rates it is reached for configurations at the mass-shedding 
limit. There are several factors which could affect the lower 
limit for rapidly rotating strange stars, such as different 
instabilities or the existence of a solid crust. The r-mode 



instability has the potential to rule out strange stars with 
small and moderate masses in LMXBs (assuming that the 
observed 1.07 kHz QPO in 4U 1820-30 is close to the 
Keplerian frequency at the innermost stable orbit of an 
accretion disk around a strange star). If the r-modc in- 
stability is not effective, then the /-mode instability and 
the viscosity-driven instability may become relevant. On 
the other hand, the dynamical instability only marginally 
affects our computed lower limits on the maximum orbital 
frequency. Our results apply to all MIT bag models. 
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